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Collective quantization of axially symmetric 
gravitating B = 2 skyrmion 
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In this paper we perform collective quantization of an axially symmetric skyrmion with baryon 
number two. The rotational and isorotational modes are quantized to obtain the static properties 
of a deuteron and other dibaryonic objects such as masses, charge densities, magnetic moments. We 
discuss how the gravity affects to those observables. 
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I. INTRODUCTION 

The Skyrme model [l| is considered as an unified the- 
ory of hadrons by incorporating baryons as topological 
solitons of pion fields, called skyrmions. The topological 
charge is identified as the baryon number B. Performing 
collective quantization for a B = 1 skyrmion, one can 
obtain proton and neutron states within 30% error Q ■ 

Correspondingly multi-skyrmion solutions are ex- 
pected to represent nuclei [1, 0, [f|. The static proper- 
ties of a B — 2 skyrmion such as mass spectra, mean 
charge radius, baryon-number density, magnetic mo- 
ment, quadrupole moment and transition moment were 
studied in detail by Braaten and Carson upon collective 
zero mode quantization Q. The results confirmed that 
the quantized B = 2 skyrmion can be interpreted as a 
deuteron. 

The Einstein-Skyrme (ES) model in which the Skyrme 
fields coupled to gravity was first considered by Luckock 
and Moss [?J ■ They obtained spherically symmetric black 
hole solutions with Skyrme hair. It is the first discovered 
counter example to the no-hair conjecture. Axially sym- 
metric regular and black hole skyrmion solutions with 
B = 2 were constructed in Q. Subsequently the model 
was extended to the SU (3) and higher baryon number 
with discrete symmetries to study gravitating skyrmion 
solutions @. 

In the Einstein-Skyrme theory, the Planck mass is re- 
lated to the pion decay c onstant f n and coupling constant 
a by M p i — / 7r y / 47r/a. To realize the realistic value of 
the Planck mass, the coupling constant should be ex- 
tremely small with a ~ O(10 -39 ), which makes the ef- 
fects of gravity on the skyrmion negligible. We, therefore, 
consider a as a free parameter and study the strong cou- 
pling limit to manifest the effects of gravity on skyrmion 
spectra. The property of the skyrmion solution could 
be drastically changed for large values of a. Certainly 



studying the effects of gravity on soliton spectra is inter- 
esting itself. Nevertheless we further attempt to give an 
interpretation to the solution as a gravitating deuteron 
or dibaryonic object and study their mass spectra and 
other static observables. 

The first step towards the study of gravitational ef- 
fects on the quantum spectra of skyrmions was taken in 
Ref. [1(3] by performing collective quantization of a B = 1 
gravitating skyrmion. It was shown there that the quali- 
tative change in the mass difference, mean charge radius 
and charge densities under the strong gravitational influ- 
ence confirms the attractive feature of gravity while the 
reduction of the axial coupling and transition moments 
by the strong gravity indicates the gravitational effects 
as a stabilizer of baryons. 

In this paper we extend the work jio} to the B = 2 
axially symmetric gravitating skyrmion and estimate the 
static properties of a deuteron or other dibaryonic ob- 
jects. We observe how strong gravity affects the baryon 
observables e.g., mass spectra, mean charge radius, 
baryon-number density, magnetic moment, quadrupole 
moment and transition moment. 

Although the Skyrme model describes nucleons with 
about 30% error, the possibility that it may provide qual- 
itatively correct description of the interaction of baryons 
with gravity can not be excluded. It is expected that in 
the early universe or equivalent high energy experiments, 
the gravitational interaction with baryons is not negligi- 
ble. We hope that our work could provide insight into 
the observations in such situations. 



II. CLASSICAL GRAVITATING B = 2 
SKYRMIONS 
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In this section, we discuss B = 2 classical regular solu- 
tions in the Einstein-Skyrme system (ES). The Skyrme 
Lagrangian coupled to gravity is defined by the La- 
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grangian 



C = Co + £s, 
lonG 



(1) 
(2) 



x [U^dpU^-^U]), 



(3) 



where is pion decay constant and e is a dimensionless 
free parameter. U describes the SU(2) chiral fields. We 
impose the axially symmetric ansatz on the chiral fields 
as a possible candidate for the B = 2 minimal energy 
configuration [|| 



U = cos F(r, 9) + ir ■ riR sin F(r, i 



(4) 



with 



= (sin9(r, 9)cosn<p, sin9(r, 6) sinrup, cos 9(r, 9)), (5) 

where n denotes the winding number of solitons, which 
is equivalent to the baryon number B. Since we are in- 
terested in B = 2 skyrmions, we consider only n = 2. 
Correspondingly, the following axially symmetric ansatz 
is imposed on the metric 



ds 2 = -fdt 2 



f 



(dr 



- 2 d9 2 ) 



r 2 sin 2 9d(p 2 (6) 



where the metric functions / , m and I are the function 
of coordinates r and 9. This metric is symmetric with 
respect to the z-axis (9 = 0). Substituting these ansatz 
to the Lagrangian one obtains the following static 
energy density for the chiral fields 

e(x,9) = ( {dxF) 2 + (^efsin 2 F) 

+ (d e F) 2 + (deQ) 2 sin 2 F + sin 2 F sin 2 9} 

I sm 9 > 



+ 



%/Isin 9 



f 



I sin 9 

(d x Fdg& - d g Fd x Q) 2 sin 2 F 



n 2 f 
I sin 2 9 



yn 2 Fsm 2 9{((9 x F) 2 
((^9) 2 + ^(9 e 9) 2 ) sin 2 F} 



;(d e F) 2 



(7) 



where dimensionless variable x = ef^r is introduced. 
The static (classical) energy is thus given by 



M = 2n^- J dxd9e(x, 9) . 



(8) 



The covariant topological current is defined by 



B> 



24tt 2 



-9 



-MU-'VvUU-'VpUUU-'VvU), (9) 



whose zeroth component corresponds to the baryon num- 
ber density 



B 



1 



-9 



■.sin 2 Fsme(d x Fd e Q - d g Fd x Q). (10) 



For the solutions to be regular at the origin x = 
and to be asymptotically flat at infinity, the following 
boundary conditions must be imposed 

d x f(0, 9) = d x m(Q, 9) = d x l(0, 9) = 0, (11) 
/(oo,0) =m(oo,0) =Z(oo,0) = 1. (12) 

For the configuration to be axially symmetric, the follow- 
ing boundary conditions must be imposed at 9 = and 
tt/2 

d e f(x, 0) = d e m(x, 0) = d e l(x, 0) = 0, (13) 
d„ f(x, ~) = d e m(x, |) = d e l(x, ~) = 0. (14) 

For the profile functions, the boundary conditions at the 
x = 0, 00 are given by 



F(O,0) = 7T, F(oo,0) = 0, 
5,9(0,0) =d x G(oo,9) = 0. 



At = and tt/2, 



d e F(x,0) =d e F(x,~) 



= 0, 

9(ar,0)=0, 9(x,|) = |. 



(15) 
(16) 



(17) 
(18) 



Since the baryon number is defined by the spatial integral 
of the zeroth component of the baryon current, we have 



B= / d 



: — (2F-sin2F) cos 9 
2tt 



(19) 



The inner and outer boundary conditions (Fq,Qq) = 
(tt,0) and (Fi,9i) = (0,tt) yield B = 2. 

By taking a variation of the static energy with re- 
spect to F and 9, one obtains the equations of motion 
for the profile functions. The field equations for the met- 
ric functions / , m and I are derived from the Einstein 
equations. We shall show their explicit form in Appendix 

m 

The effective coupling constant of the Einstcin-Skyrme 
system is given by 



a = AnGf 2 



(20) 



which is the only free parameter. 

We use the relaxation method to solve these nonlinear 
equations with the typical grid size 100 x 30. In FigfTJ 
the profile functions for a = 0,0.04,0.08,0.126 is pre- 
sented. No solution exist for a > 0.127. Also, the metric 
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functions are shown in Figs[2]-|4] In Fig[5j we display the 
metric functions at 9 — 0, 7r/4, tt/2 as a function of ra- 
dial coordinate. Fig[6] shows a dependence of the static 
energy M in unit of /w/e. In Fig. [7] the energy densi- 
ties defined in Eq.([7]) are plotted. The baryon densities 
b = ^—gB° are shown in Fig. [5) 



III. COLLECTIVE QUANTIZATION SCHEME 
FOR THE B = 2 SOLITON SOLUTION 



In this section, let us try to give an interpretation to 
the B = 2 skyrmion solution obtained in the last section 
as a baryonic object by assigning quantum number of 
spin and isospin to it. We employ the semi-classical zero- 
mode quantization following Ref. [||. Let us introduce 
dynamically rotated chiral fields 

U(r,t) = A{t)U{r')AHt), r' = R(B(t))r, (21) 



and 



Rij (B) 



TrfaBTjB*), 



(22) 



where A(t) and B(t) are the time dependent SU(2) ma- 
trices generating the isospin and the spatial rotations. 
Substituting ([21) into the Lagrangian |3]), one finds 



L = -M + -a i U lj a j - <i,\Y,,h, + -biVijbj, (23) 

where M is the static energy of the B = 2 skyrmion 
given in Eq.©. The Lagrangian is quadratic in time 
derivatives 



iYtt,A^A , bi = iTtt,BB\ 



(24) 



The moments of inertia tensor Uij , Vij and are ex- 
pressed in terms of the chiral field U as 



U, 



,, — J ■ ■ L ^x 2 sin0{Tr(C/t[^,C/][/t[^,[/]) 

+ g kl Tr[U^d k U, U]] [U^d t U, pt[S [/]] }] ,(25) 

W lj = U l] {[^,U] -> i{x'times^) 3 U)} , (26) 
V ij =W ij {[f ,U] -> i(x'timesV) z U)} . (27) 

Body-fixed isospin operator Ki and angular momentum 
operator Li are defined as a canonically conjugate to 
and bi 



dL_ 

dL 
db~ 



Uijdj 



Wtjbj, 



Vijbj. 



(28) 



These operators are related to the usual coordinate-fixed 
isospin Ii and spin Ji by the orthogonal transformation, 

Ii = -Rij {A)Kj , Ji = -Rij (B) T Lj. (29) 

The commutation relations for these operators are 

[Ki.Kj] — ieijkKk , [Li,Lj] = ieijkLk, 

= i^ijkh , [Ji,Jj] = itijkJk- (30) 

These operators satisfy the Casimir invariance by using 
Eq.~ 



K 2 = I 2 , L 2 = J 2 . 



(31) 



The symmetry of the classical soliton induces the follow- 
ing conditions for the inertia tensors 

U u = U 2 2 , V u = V22 , W n = W22 = 0, 

^33 = 4£/ 33 , W 33 = 2U 33 . (32) 

Thus it is sufficient to calculate only Uu , U 33 and Vu- 
Inserting Q into Eqs.fJ^S]) and (|2T|) one obtains 



Uu 



77^ 



dxdO 



4/ 



^x 2 sin 9 sin 2 F(l + cos 2 0) 



U 



33 



he 



+ ^ sin 6» sin 2 F{{1 + cos 2 Q){x 2 F 2 x + F 2 g ) 
+ sin 2 F cos 2 Q(x 2 2 x + Q 2 e ) 
+ 7 ^sin 2 Fsin 2 e}j, 

dxd9 \^x 2 sin 9 sin 2 F sin 2 6 
2& sin 9 sin 2 F sin 2 6{a; 2 F 2 + F 2 g 
+ sin 2 F(x 2 Q 2 x + <d 2 g)}], 



(33) 



(34) 



Vu = 



77? 



dxd9 



^x 2 sin (9(F 2 + B 2 g sin 2 F 



-n 2 cot 2 6» sin 2 F sin 2 6) 



+ v^2 sin sin 2 F{(F x 6, e - F fi Q. x ) 2 



f 



1 2 

,2/77.2 



2 (F 2 + 6 2 x sin 2 F) cot 2 9 sin 2 9} 
+ SI(cos 2 e+f) 

x (F 2 e + e 2 e sin 2 F) sin 2 F sin 2 el . (35) 



From Eqs. (j28]) and (1321) we derive the constraint 

(2F 3 + L 3 ) \phys) = 0. (36) 



Inserting the body- fixed operator (|28p , Casimir invariant 
((31]) and constraint f|36|) into the Lagrangian ([23]) one can 
get the Hamiltonian operator as 



H = M 



J 



2Uu 2Vu 

1 1 4 

f/33~ _ t/n ~ t^T 



K 2 

If- (37) 
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a = 0.000 F(P,Z) 3 a = 0.000 @(P,Z) 1.4 






FIG. 1: The profile function F and 9 in the cylindrical coordinate with a = 0.000,0.040,0.080,0.126 are shown. We use 
dimensionless variables P = ef n p and Z — ef n z. There exists no solution for a > 0.127. 



5 



a = 0.040 f(l,Q 




a = 0.040 m(i,Q 



0.95 
> 0.9 
^0.85 h 
5 0.8 

0.75 



m 
- ii 










B V 

o.6 V r~ 






0.4 V 

0.2V 



* — trr 1 °^ ! 


te '' VA ''Of 



a = 0.080 m (%Q 



C.995 
0.98 

0.965 
0.95 



<^0.85 
0.8 
0.75 
0.7 



a = 0.126 /fcO 




FIG. 2: The metric function / in the cylindrical coordinate 
with a = 0.040, 0.080, 0.126 is shown. We use £ = P/(l + P) 
and £ = 2/(1 + Z), instead of dimensionless variables P = 
e/,rp and Z = ef^z. 



The corresponding energy (mass) eigenvalues are 
»(*+!) , J'O' + I) 



E = M 



2Un 
1 

033 



2Vii 
1 4 



T 



(38) 



where i(i + 1), j(J + 1) and k are the eigenvalues of the 
Casimir operators (|3"Tj) . and the operator .K3, respec- 




FIG. 3: Same as Fig[2]for the metric function m. 



tively. The parity is defined by the eigenstate of the 
following operator 



P=e 



(39) 



which means that the parity is (+) for even k and (— ) 
for odd k. 

Taking into account the Finkelstein-Rubinstein con- 
straints for the axially symmetric soliton [6], one finds 
that some combinations of are not acceptable. The 




FIG. 4: Same as Fig[2]for the metric function 

allowed states are 0, [HJ 

IM3O) \jj3O) , provided i + j is odd, 
(for k = Q) 

^[|"'3*> \J3S ' 2 «> - l«3 - «) Iij32/t>]. 

(tarn =!,■■■ ,min{z,[j/2]}) (40) 

Therefore k > is allowed only if i > 1 and j > 2. 
We show quantized energy spectra for various values 



FIG. 5: The metric function /, m and / for — 
0, 7r/4, 7r/2 as a function of radial coordinate with a = 
0.040,0.080,0.126 are shown. We use X = as/(l + »), 
instead of dimensionless variable x = e/»r. 



of a in Table HI The mass difference from the classical 
energy is shown as a function of a in Fig. [5] In Table U 
(and FigO, we show the results for i,j < 3 with k = 0, 
and also some excited states for k < 2. We present a spec- 
troscopic classification 2s+1 Lj in Table Q] For n = 0, the 
system is in S- state, and then j equals to the intrinsic 
spin s. For k > state, the orbital angular momentum 
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^ 55 



0.02 0.04 0.06 0.08 0.1 0.12 



FIG. 6: The coupling constant dependence of the classical 
B — 2 skyrmion dimensionless mass derived from Eq.©. 



is chosen to be its lowest value so that it is consistent with 
the quantum number j , the parity and the fact that s < 3 
in the six-quark picture. It is seen that the mass differ- 
ence increases monotonically with increasing a. Fig. [1] 
(and also Figs l7l8p implies that the strong gravity makes 
the size of the skyrmion smaller which makes the inertial 
moment smaller, resulting in increase in the mass differ- 
ence. In the collective quantization, the skyrmion can be 
quantized as a slowly rotating rigid body and the mass 
difference is interpreted as a consequence of the rotational 
kinetic energy. Thus the gravity works for increasing the 
kinetic energy of the skyrmion. In the naive quark model 
picture, the mass difference is ascribed to the hypcrfinc 
splittings. The increase in the mass difference may imply 
that due to the reduction of the distance between quarks, 
the effects of the hyperfine splittings become dominant 
by the gravity [HI]. 

From the data of the Table[l] one can straightforwardly 
obtain the mass value in MeV unit, by multiplying f„/e 
for any parameter choice (/ w ,e). With f„ = 108 MeV, 
e = 4.84, which is used in Ref. [6], our obtained mass 
value in the case of a = is approximately 6% lower 
than the value obtained in Ref. @. There are two pos- 
sible reasons for that. Firstly, the Lagrangian in Ref. @ 
includes the mass term and hence it slightly enhances the 
mass value. Thus we also performed the same analysis 
including the mass term, which reduced the error from 
6% to 0.5% in the mass. Secondly, the coordinate system 
they used is cylindrical while wc adopted the spherical co- 
ordinates. Since the cylindrical coordinate requires much 
larger number of grid point to obtain torus-shape solu- 
tions than the spherical, we suspect that their results are 
not fully convergent, producing a slightly larger mass. 
IV. ELECTROMAGNETIC PROPERTIES 

In this section we shall investigate the gravity effects 
to the various observables, i.e. the mean charge radius 



•2\ 



1/2 



, magnetic moment /id, the quadrupole moment 



with the process 7c? — > 1 So- 

Let us introduce the electromagnetic current in the 
Skyrme model jjf m \x) which consists of the isoscalar 
and isovector part, 

J^\x)= 1 -B tl {x)+Il{x), (41) 

where B^(x) is baryon current density given in Eq.Q 
and I 3 (x) is the third component of the isospin current 
density, defined by 

xp^d^U-^U]). (42) 

Inserting the dynamical field U in Eq. (j2Tj) into Eq. ([4Tj) 

gives electromagnetic current operator Jp Cm ' . 

To estimate the expectation value of these quantum 
operators, we describe the quantum spin states (140|) in 
terms of the products of the rotation matrices 



(A\n 3 k 3 ) 
(B\jj 3 h) 



2i + l 
2j + l\* 



D\it 2 A%, 



(43) 
(44) 



where D l (A) mt7n i is well known Wigner D function. The 
computation can be easily performed by the following 
integration formula 14 1 



J dB D\B)l im ,D3(B) m2m ,D k (B) 



2k -\- 1 im i3 m 2 im' 1 jm' 2 '' 



(45) 



where C^ i 3 jm2 is a Clebsch-Gordan coefficient. 

1 /2 

The deuteron charge radius (r 2 ) d is defined as the 
square root of 



(r 2 ) d ^(d\l d 3 rr*J^\r,t)\d) 



(46) 



where | d) represents spin state of the deuteron. In this 
case, as only the isoscalar part of Jq CI1 ^ contributes to 
the matrix element, the integration is straightforward. 

{r2)d = Vt? I dxd9 ^~9 x2 B (r,9) 



dxdOx sin F sin 6 



x(d x Fd e e-d e Fd x Q). (47) 



Q, and the transition moment fid^np, which is associated 



FigfTOlshows the a dependence of the mean charge radius. 
Due to the attractive effect of the gravity, it decrease with 
increasing a. 
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TABLE I: The dimensionless quantized energy spectra up to i,j < 3 for k = 0, and also some excited states up to k < 2. In 
the particle classification, we give the baryonic descriptions if available, otherwise we only show its spectroscopic classification 

2s + l r 



Classification 


i 


j 


K 


Parity 


a = 


a = 0.040 


a = 0.080 


a = 0.126 


classical skyrmion 










69.7195 


63.6722 


56.9827 


44.4862 


Deuteron ( a Si) 





1 





+ 


69.7227 


63.6758 


56.9866 


44.4912 


NN (^o) 


I 








+ 


69.7241 


63.6774 


56.9886 


44.4939 


rp 2 ) 


I 


2 


1 




69.7294 


63.6828 


56.9943 


44.5006 


NA ( b S 2 ) 


1 


2 





+ 


69.7339 


63.6881 


57.0003 


44.5089 


NA (' 5 S 2 ) 


2 


1 





+ 


69.7368 


63.6913 


57.0041 


44.5143 


( J P 2 ) 


2 


2 


1 




69.7387 


63.6932 


57.0059 


44.5160 


AA ('Ss) 





3 





+ 


69.7391 


63.6935 


57.0062 


44.5162 




I 


3 


1 




69.7392 


63.6935 


57.0060 


44.5156 


AA OSo) 


3 








+ 


69.7475 


63.7032 


57.0176 


44.5324 


(°D 4 ) 


2 


4 


2 


+ 


69.7479 


63.7024 


57.0153 


44.5262 


( ;, Pa) 


2 


3 


1 




69.7485 


63.7039 


57.0177 


44.5310 
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= 0.000 Baryon density b(P,Z) 



a = 0.040 Baryon density b(P,Z) 




FIG. 8: 

The baryon density b = y/—gB° in the cylindrical coordinate with a = 0.000,0.040,0.080,0.126 are shown, in terms of the 
dimensionless variables P = ef^p and Z = ef n z. 



The isoscalar part of the magnetic moments is ex- 
pressed in term of electromagnetic current as 

Ai= \J d*r^e im j£ m \r,t). (48) 

Inserting the dynamical field U in Eq. ([2T1) into Eq. (|48|) . 
one finds the magnetic momentum operator as 

Ai| 7=0 = \{RiA B ) T ^jkak +n' jk b k }, (49) 

where 

Mjfc = j-j y J d 3 x^Eji m xi^Tr[U^[^T k , U]C m ], (50) 

1 f 1 1 

/4* = T~Fv2 / d 3 x-£ji m xi-Tr[W(-ix x V) fc E/ 

x<7 m (as)],(51) 

and C m is 

C ro (x) = -^e mnp (U^d n U){rfd p U). (52) 



The i?y is the rotation matrix ([22]) and di and &j are 
defined as Eq.fJM]). The anti-commutator relation in 
Eq. (|48p guarantees fa to be Hermitian operator, fiij 
and n'j k are diagonal and furthermore they satisfy fin — 
M22 = , /i' n = p! 22 and ^ 33 = — 2/133. With these 
relations and using the body-fixed operator (f2"B"|) , the 
coordinate-fixed operator (f2"9")l , the relation of the mo- 
ment of inertia components 1|32[) , and the constraint of 
Eq.([3BD in Eq.([lg|l. one can get 

^Ii— = — V~"^ 1 ^ erms P ro P or tional to -K3. (53) 
Therefore we only use the component of // 1:L , 
Mn = - 7 , 2 / dxdBj~gx 2 (cos 2 + 1)Bq(x, 9). (54) 

And the Vn is moment of inertia given in Eq. (|35|) . The 
moments fid of deuteron is defined as the expectation 
value of the with J3 = 1 

W = (h) = (55) 
1 
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a 

FIG. 9: The coupling constant dependence of the dimension- 
less mass difference from the classical energy are shown. For 
a > 0.127, there exists no solution. We give the baryonic 
descriptions if available, otherwise we only show its spectro- 
scopic classification ( 2s+1 Lj)i, whose subscript i means eigen- 
value of the isospin operator (I31|l . 
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1.2 ' 1 ' 1 1 ' 1 — 

0.02 0.04 0.06 0.08 0.1 0.12 
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FIG. 10: The coupling constant dependence of the dimen- 
sionless mean charge radius (|47)l . 
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FIG. 11: The coupling constant dependence of the dimension- 
less magnetic moment (|55|l . 



In FigQj] a dependence of the [id is shown. It decreases 
with increasing a, but the effect of gravity is more evident 
compared with that of B = 1 10' . In the non-relativistic 
nuclear physics point of view, the magnetic moment of 
the deuteron consists of the sum of the magnetic moment 
of the proton plus neutron and some correction terms re- 
lated to the D-state probability of the deuteron. This 
result indicates that the assumption that the deuteron is 
almost S-state (contributions to the other states are only 
~ 5%) may no longer valid under the strong gravitational 
field. Therefore, we speculate that the effect of gravity is 
apparent on the change of such non-central components 
of the nuclear force. This will be more evident by exam- 
ining the quadrupole moment. The quadrupole moment 
is given by 

Qa = J <l l r, —.:>„;,■, - r 2 5ij) J (em) (r, 9), (56) 

which is the same as the magnetic moment. Inserting 
(|2ip into Eq. (|56[) one obtains the quadrupole momentum 
operator 

Qij\ I=0 = Ria(B) T Q a bRbj(B), (57) 

where 

Qab = , \ , 2 ( d 3 xy^g(3x a x b - x 2 5 ab )B a (x, 9). (58) 

Qij satisfy Qn = Q22- And the symmetry relation for 
the quadrupole moments reduces the expression to 

Qij\i=o = Q^l^(B) T R 3j (B) - itfy]. (59) 
Thus we only need the component of Q 33 

Qa3 = -r-T-^ I dxde^/^jx 2 (3 cos 2 9- l)B°(x,6)(ffl) 
{eUY J 
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FIG. 12: The coupling constant dependence of the dimen- FIG. 13: The coupling constant dependence of the dimen- 
sionless quadrupole moment sionless transition moment ([551). 



The deuteron quadrupole moment is defined as the ex- 
pectation value of Q33 with j 3 = 1 

<2=(<933>, 



^(d,j' 3 = l\B% 3 (B)R 33 (B)\d,j 3 = l) - ~ 



(61) 

Now let us evaluate the matrix element. The rotating 
matrix R 33 is represented as R 33 (B) — D 1 (B) Q0 with 
the Wigner D function. The product of two £>-functions 
can be expanded in the following series fl4j 



D Jl (B)m 1 n 1 D J2 {B)m 2 n 2 



J1+J2 



MNCj^ lihN2 . (62) 



- Y C JiMiJ2M 2 D ' J ( B ) 
J=|Ji-Ja| MN 

The matrix element of D 1 (B)qo derived from Eqs. (|44l 
and (|45[) and the relation of the D l (iT2)om = (— 1) $ 



m.Oi 



{j'j' 3 1' 3 \d\B) Q0 \jj 3 1 3 ) = -(^^)*C^ 10 CJ* 10 (63) 

The deuteron state is the one with j 3 = 1. Therefore 
inserting i = 0, j = X, k = and using the constraint in 
Eq. (j36)) . the matrix element in Eq (|6"l"j) is evaluated as 



Q = -IQs3- 

5 



(64) 



In Fig[T2] a dependence of the Q is shown. As in the 
case of the magnetic moment, the quadrupole moment 
significantly decreases with increasing a, which suggests 
that the D-state probability is strongly affected by the 
gravity. 

The transition moment fid—>np f° r photodisintcgration 
of the deuteron into the isovector 1 So state is defined 
by the magnetic moment operator jx 3 . It is evaluated in 
terms of matrix element between the j 3 = state and 
i 3 = state. 



^Sq,^ = 0\fi 3 \d, j 3 = 0). 



(65) 



Inserting the dynamical held in Eq. (PiTj) into the isovector 
part of fii, one can obtain 



1 



£i| J=1 = -^R3j(A)W jk R ki {B) 



(66) 



where Wj k is the moment of inertia tensor defined in 
Eq. tpp]) . From Eq. (f52]) . one can see that W33 is the only 
non-zero component and W33 = 2U 33 , and hence 



/'3 



7=0 



-U 33 R 33 (A)R 33 (B) 



(67) 



Evaluating the matrix elements of (f65|) using (163]) and the 
element of the isospace which is derived in the similar way 
in Eq. (|6"5]) . one can get 



(68) 



In Fig |13l the coupling constant dependence of fj,d-* n p is 
shown. As expected, it also decreases with increasing a. 
Since the strong gravity reduces the transition moment 
significantly, it may be possible to determine the grav- 
itational constant by observing the variation in fj,d~*np- 
It is interesting that the decay rates are reduced by the 
gravitational effects whether the interaction is strong or 
electromagnetic, which means the gravity works as a sta- 
bilizer of baryons. 



V. CONCLUSIONS 

In this article we have investigated axially symmetric 
B = 2 skyrmions coupled to gravity. Performing col- 
lective quantization, we obtained the static observables 
of the dibaryons. The rotational and isorotational modes 
were quantized in the same manner as the skyrmion with- 
out gravity. It was shown how the static properties of 
dibaryons such as masses, charge densities, magnetic mo- 
ments were modified by the gravitational interaction. 

The dependence of the energy density and mean radius 
on the coupling constant showed that the soliton shrinks 
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as a become larger which reflects the attractive feature 
of gravity. The mass difference between dibaryons also 
becomes larger as increasing a. These observations can 
be interpreted that shrinking the skyrmion reduces the 
inertial momenta and hence induces the large mass dif- 
ference between dibaryon spectra. 

In the collective quantization, the skyrmion can be 
quantized as a slowly rotating rigid body and the various 
dibaryon spectra are regarded as the rotational bands 
of the classical skyrmion. Thus the large mass difference 
means that the gravity works for increasing the rotational 
kinetic energy of the skyrmion. The magnetic moment 
is reduced significantly as increasing a compared to that 
of B = 1. We also calculated the quadrupole moment. 
For the transition moments of the deuteron, we found 
that the gravity works as a stabilizer. Thus, it may be 
possible to determine the gravitational constant by the 
measurement of the various decay rates of the dibaryonic 
objects. 

Throughout the paper, we consider a as a free param- 
eter and studied skyrmion spectra in the strong coupling 
limit. We found that the effects of gravity on the ob- 
servables estimated here are manifested only in such a 
large coupling constant. Some theories such as scalar- 
tensor gravity theory fill ] and theories with extra dimen- 
sions discuss the time variation of the gravitational con- 
stant 16]. There may have been an epoch in the early 
universe where the gravitational effects on nucleons were 
significant. 

Let us note that in Ref.[I3| the authors computed the 
quantum correction to the B = 2 skyrmion by con- 
sidering the moduli space of instanton-generated two- 
skyrmions in the attractive channel which has a larger 
dimension (.Mio) than the moduli space of the deuteron 
(.Ms), which greatly improved the deuteron observables. 
Thus, going beyond the collective coordinate approxima- 
tion provides a much greater correction to the observables 
than gravitational effects. 

As a future work, the analysis of the electromagnetic 
form factors of the deuteron is now under consideration. 
The SU(3) extension and the analysis of dihyperon cou- 
pled to gravity will be also interesting. 
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APPENDIX A: FIELD EQATIONS 

For the numerical calculation, we employ the radial 
coordinate 



1 + x 



(Al) 



to map infinity to \ — 1- The regularity requires the 
condition 



m( X ,0)=J(x,0). 



(A2) 



Following the work of Ref. [9( we introduce the metric 
function g as 



g(x,o) 



i(x,o) ' 



(A3) 



Then the boundary conditions for g are given by 



g(0,6) = l , g(oc, 



KX,0) = 1 , g( x ,-)=Q. 



(A4) 
(A5) 



The filed equations for the profile function F(x,0) and 
0(x, 0) are derived as 



SF J /| VI sin 
X 2 (l - X ) 2 (F XX + l -$F x ) + 2 X (1 - X?F X 



+F : ee +F,eC-£ + cot( 



-{X 2 (l " X) 2 © 2 X + e 2 e }sinFcosF 
-7^^sinFcosFsin 2 G 



4^ sin 2 F 



'la 
- / 



(i-x) 4 {(^ + ^-^)e, 9 

+ cotF{F )X Q t0 - F e Q >x )}{F iX e, e - F e 6, x ) 
+{e, e (F, xx Q, e + F, x e, x6 - F, x gQ, x - F, e ®, xx ) 
-Q !X (F <x gQ >e + F tX Q fi e - i>9, x - F fi Q, x e)} 
x (1 - x f ~ 2(1 - x ?®fi{F, x Q,e - F, e Q, x ) 



in 2 f ■ 2 n ■ 2a 

- , . i a sin b sin H 

/ sir 8 



(i-x) 4 {^ 



+F 2 X cot F + 2F x 6 iX cot 6 + (if - l -§)F t 



i 



2 sin F cos FQ 2 X } - 2(1 - x?F tX 
^^{F^ee + F 2 g cot F + 2F e Q,e cot 6 

2 



(if - ^ -cot0)F„ -26^sinFcosF} , (A6) 
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'S^C S \ 1 



<5© J f 2 sin (9 sin 2 F 

= x 2 (i - x) 2 (e xx + % e, x ) + 2 X (i - x) 2 e, x 
+e, e0 + e, e (^ + cote) 
+2{ x 2 (i - x ) 2 F x e, x + F, s e, fl } cot f 

-OTsinOcose 

+4^[(i-x) 4 {(^ + ^f-^r + cot ^x 
)F e }(F x e, 9 -F e e, x ) 



2/ 



+{F x (F xe e, e + F x e, ee - *>e iX - F fi e iXe ) 
-F e (F xx e, e + F x e, xe - F x ,e, x - F 9 e, xx )} 

x (1 - x) 4 + 2(1 - X fF.g{F x Q.g - F fi Q tX ) 

(i-x) 4 {e,, 



. sin 2 Fsin 2 

/ sin z # 



.XX 



+e 2 x cot e + 4F x e, x cot f 
+(^f-^)0,x}-2(i-x) 3 e, x 
+ (i^xl!{e ee + ®z e C ot e + 4F e e, cot f 



(lz.x) 

X 

+ (^-^-cot0)e, 9 } 



i^L sin 6 cos 6 



(1 - X) 4 F 2 X + ^F 2 



• (A7) 



From Einstein equations, the field equations for the met- 
ric functions are derived. We diagonalized those equa- 
tions to the 2nd-derivative of x and for each metric 



function 



j|x 2 Goo + x 2 G\\ + G 22 + 1 J^i g G33 
-2a[fr 2 T m + f 2 T n + T 22 + j^gT 33 

= X 2 (l - X) 2£ T X 2 (l - X) 2 (^f) 2 + 2 X (1 - X) 2 ^f 

+if- - (^) 2 + cote^ 
+i Y 2 (i - y) 2 ^ 1 -^- + ±L°-h!- 

t 2 X v 1 X) f l ^ 2 / l 
-a [2^(1 - x) 4 (F x e, 9 - F e 9 x ) 2 sin 2 F 
+ ll^{(l-x) 4 (^ x + e 2 x sin 2 F) 
1 ii =^(F 2 + 2 g sin 2 F)} sin 2 F sin 2 e 



4P 



(A8) 



I sin 2 9 



2a 



m rp 
I sin 2 1 33 



X 2 (l-X) 2 ^ + ^X(1-X)(l-2X)^ 



2 A V x A7 V m / 2 m 2 V m ) 

+ j[x 2 (l-x) 2 (F 2 x + e 2 x sin 2 F) 



-F 



<0 t Q 2 g sin 2 F - J f^ J sin 2 F sin 2 9 



-4^(1 - X ) 4 (F x e, e -F 9 e, x ) 2 sin 2 F 



4n J / „;„2 
I sin 2 

(i-x) 



sin 2 Fsin 2 6{(1 - x) 4 (^, x + © 2 X sin 2 F) 



+ ^(F 2 +e 2 ,sin 2 F)}] 



a; 2 Gii + G 22 - 2a 



(A9) 



x 2 T n + T 22 

= k 2 (l-x) 2 ¥-3X 2 (l-x) 2 (¥) 2 

+X(1 - X)(| - xY-f + \C-ff + cot 8 l -f 

-a[2i(l- X ) 4 (F x e, e -F e e, x ) 2 sin 2 F 



-~/^sin 2 Fsin 2 el . 

2/ sin' e J 



(A10) 
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